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Abstract 

The coupling of electronic and vibrational motion is studied by two canonical transforma- 
tions namely normal coordinate transformation and momentum transformation on molecular 
Hamiltonian. It is shown that by these transformations we can pass from crude approxi- 
mation to adiabatic approximation and then to non-adiabatic (diabatic) Hamiltonian. This 
leads to renormalized fermions and renormalized diabatic phonons. Simple calculations on 
H2, HD, and D2 systems are performed and compared with previous approaches. Problem 
of electronic quasi- degeneracy is discussed. 
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I. INTRODUCTION 



Many atomic systems (e.g. molecules and/or crystals) are systems with interacting 
electrons and nuclei and can be thus described by Schrodinger equation 

= Em (1) 

In general case the number of degrees of freedom of such systems with coulomb interaction is 
too large and even using high-performance computers to solve this equation becomes impos- 
sible. Therefore the only possibility to solve eq. @j is to suggest some approximations 
The most important approximation and also most often used is the Born-Oppenheimer (BO) 
and the adiabatic approximation. This approximation is based on the fact that masses 
of nuclei are 1870 times heavier than the masses of electrons. This leads to the idea of 
potential energy surface. Beside the many attemps to go beyond the BO approximation and 
many different approaches cartain facts are not clear completely. In this paper we decided 
to study the coupling of electronic and vibrational motions by two canonical transforma- 
tions namely normal coordinate and momentum transformations. Our approach is similar 
to quasiparticle transformations often done in solid state physics. In order to make our 
approach more transparent we repeat here main features of adiabatic approximation. We 
follow the arguments of recent Kutzelnigg's paper 
Let us write the total molecular Hamiltonian as 

H =TN{R) + ENN{R) + HEN{r,R) + HEE{r) (2) 

where Tn{R) is the kinetic energy of nuclei, E^j^{R) is the interaction between nuclei, R 
denotes nuclear coordinates, r denotes electronic coordinates, and HEN{r,R) + HEE{r) is 
the electronic Hamiltonian 

HENir, R) + HEEir) = TE{r) + UEN{r, R) + HEE{r) (3) 

where Tjv(r) is the kinetic energy of electrons, UEN{f, R) is the electron-nuclei interaction 
term and HEE{f) represents the electron-electron interaction. In Born and Huang 
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approach the total wavefunction depending on the nuclear coordinates R and the electronic 
coordinates r is expanded as 



where ipk are a complete set of known functions of r that depend parametrically on 
the nuclear coordinates R and where the Xk{R) are regarded as unknown. The ipk are 
conventionally chosen as a set of eigenfunctions of the clamped nuclei (CN) Hamiltonian, 
but this is not necessarily a good choice, since this is usually not complete (without, the 
corresponding continuum functions, which one cannot include anyway). Both the BO and 
adiabatic approximation can be based on choosing a single term in eq. @ 



This is referred as BO ansatz. This ansatz is taken as a variational trial function. Terms 
beyond the leading order in m/M are neglected ( m is the electronic and M is nuclear mass, 
respectively). The problem with expansion (1) is that functions ?/'(r, R) contain except bound 
states also continuum function since it includes the centre of mass (COM) motion. Variation 
principle does not apply to continuum states. To avoid this problem we can separate COM 
motion. The remaining Hamiltonian for the relative motion of nuclei and electrons has then 
bound state solution. But there is a problem, because this separation mixes electronic with 
vibrational coordinates and also there is a question how to define molecule-fixed coordinate 
system. This is in detail discussed by Sutcliffe |^ . In the recent paper by Kutzelnigg |^ 
this problem is also discussed and it is shown how to derive in a rigorous manner adiabatic 
corrections using so called Born-Handy ansatz. There are few important steps to arrive at 
formula for a diabatic corrections. Firstly, one separates off COM motion. Secondly, (very 
important step) one does not specify the relative coordinates (which are to some extent 
arbitrary). In this way one arrives at relative Hamiltonian Hrei [H with trial wavefunction 
"^rei- If we make BO ansatz 




(4) 



k 




(5) 




(6) 
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where pik, p^y are non-specified relative coordinates and ip is chosen as a solution of the CN 
Schrodinger equation. The adiabatic correction AE take very simple form 

AE = -]^<m\Y.M,Vl\^> (7) 



This formula was used by Handy [|I^]. It was used previously also by Sellers and Pulay 
TT| . (See also Davidov for derivation). For practical calculation the identity 



f f dip \ ^ 
can be used. 

Note that practically in any textbook the validity of BO approximation is justified only 
when 

< 1 (9) 



\En — En 



when uj is the frequency of harmonic vibrations arround the point Rq. 
The aim of this paper is twofold: 

i, We show how starting with molecular Hamiltonian in crude adiabatic representa- 
tion we arrive at adiabatic Hamiltonian by performing canonical transformations which mix 
together the electronic and vibrational motions (through normal coordinates). We derive 
simple formulae for adiabatic corrections, similar to eq. (^. 

ii, We generalize canonical transformations (through momenta) arriving at non-adiabatic 
Hamiltonian. We introduce the idea of quasiparticles (renormalized electrons) and present 
the formulae how to obtain the "orbital energies", "correlation corrections" and non- 
adiabatic frequencies for these quasiparticles. Finally, we perform some simple model cal- 
culations to demonstrate how the method works. 

II. THEORY 

Let us start with electronic Hamiltonian (Bl) which we denote 
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HEN{r, R) + HEE{r) = h + p\ (10) 

where h is the one-electron part representing the kinetic energy of the electrons and electron- 
nuclear attraction term, and f ° is the two electron part of the Hamiltonian corresponding 
to electron-electron repulsion term. For the purpose of diagrammatic many-body pertur- 
bation theory it will be efficient to work in second quantization formalism. The electronic 
Hamiltonian (ITUp has the form 



PQ 

+ - < PQ\v^\RS > a+a+asan (11) 

^ PQRS 

where cbp{aQ) is the creation (annihilation) operator for electrons in the spinorbital basis 
|P >, IQ >, ■ ■ ■ • If we apply the Wick theorem to (0) we can write this equation as 



H 



EN 



+ Hee= E^// + ^E( 



^hlj-^'ljjl) + 



J2hpQN 

PQ 



IJ 

a^aq 



E 





PIQI 



VpiIQ 1 N 



+ 9 E ^PQRSN 
^ PQRS 



PQI 

apCLQCigCi fj 



apttQ 



(12) 



where 







ABABK'^ABBA 



denotes the coulomb (exchange) integral. One possibility is to work 



the within crude representation in which the spinorbital basis |P >, |Q >, ■ ■ ■ is determined 
at some fixed (equilibrium coordinate Ro). Note that Hamiltonian (^) has 3A^ — 6 degrees 
of freedom (in fact 3N degrees of which 6 are zero). Hamiltonian (|T^ has only bound-state 
solutions. Let us divide individual terms of the Hamiltonian ([12|) into two parts. Namely 
caculated at point Ro and the terms which are shifted with respect to term at Ro (we use 
prime to denote these terms). The electronic Hamiltonian (|l^) can be rewritten as 



Hen + Hee = E'^scf + ^scf + E 



apap 



+ E h'pQN 
PQ 



ttpaq 



+ 9 51 '"PQRS 
^ PQRS 



N 



O'pCtiQO'S^R 



(13) 



where E^^jp is the Hartree-Fock energy calculated at the point Rq , and h'^^p is the shift 
in the Hartree-Fock energy with respect to the point other than Rq. The same is true for 
one-particle operator of (]T3|), where ep are the one-particle Hartree-Fock energies calculated 
at point Rq. The correlation operator is not changed because it does not depend on nuclear 
coordinates R . For the notation see ||12|. Let us perform the Taylor expansion for the 
energies -Eatat and uscf around the point Rq. 



E 



NN 



/ . ^NN 
i=0 



and 



(14) 



Uscf 



u 



(0) , ' 

SCF + Ug(jp 



E 

i=0 



u 



(i) 

SCF 



Using (jl^) and ([ISD we can rewrite our Hamiltonian ([TsD in the form 



H = EZ + 4cF + T.^^r( btbr + ^ j + uPcF + E ^pN 



CipCLp 



+ J E ^PQRSN 
PQRS 



CLpCLQttsO'R 



where 



+ -^ATAT E^ + Ugfjp 



U 



(2) 
5CF 



E^PQ^ 
PQ 



apttQ 



T J- , „(2) 



(15) 



(16) 



(17) 



ujr is the frequency of the harmonic oscilator and [b) are boson (phonon) creation (annihi- 
lation) operators. In order to use the perturbation theory we have to split the Hamiltonian 
(p^) onto the unperturbed part Hq and the perturbation H' 



H 



(18) 



Due to the crude approximation, we can partition the Hamiltonian (|T^) in the following 



way 



Ho = + 4°^^ + ^ epN \a+ap] + hur f + ^ ) (19) 



p ~ ' r 



and 

R' = H'j, + H'p + H'j. (20) 

Where H' contains all the terms in (|l^) except (19). In eq. (|l^) all quantities were defined 
through the cartesian coordinates. For further purposes it will be natural to work in normal 
coordinates {-Br}- The normal coordinate in second quantized formalism is given as 

Br = hr + h+ (21) 

If we transform Hamiltonian (ff^) into normal coordinates we arrive at the following expre- 



sions 



H = E% + 4°^^ + J2 [atap] + ^ hu., (b+K + ^ 

p ^ V z 

+ ^ E ^PQRS^ [at.a^asaR] { = H'^} 

^ PQRS 

oo [n/2] 

ji=l,(ri^2) fc=0 

oo [n/2] 

+ E E E u%-''\B(-^'^N[ataQ] { ^ H',} (22) 

n=l k=0 PQ 

The term H'^ is the electron correlation operator, the term Hp corresponds to phonon- 
phonon interaction and Hj corresponds to electron-phonon interaction. If we analyze the 
last term Hj we see that when using crude approximation this corresponds to such phonons 
that force constant in eq. (|l^) is given as a second derivative of electron-nuclei interaction 
with respect to normal coordinates. Because we used crude adiabatic approximation in which 
minimum of the energy is at the point Rq , this is also reflected by basis set used. Therefore 
this approximation does not properly describes the physical vibrations i.e. if we move the 
nuclei, electrons are distributed according to the minimum of energy at point Ro and they do 
not feel correspondingly the R dependence. The perturbation term Hj which corresponds 



to electron-phonon interaction is too large and thus perturbation theory based on sphtting 
given by eq. (|19|, ^) will not converge . Natural way to improve this situation will be to 
use basis set which is generally R dependent. We can do this in second-quantized formalism 
in a way that we pass from electron creation (annihilation) operators ap (ag) which act on 
Ro dependent basis set to a new fermion creation (annihilation) operators ap ( aq) which act 
on R dependent basis. Similar transformation was studied for solid state theory by Wagner 



1^ , who also discuss the convergency properties of adiabatic approximation [|15|. This we 
can achieve by canonical transformation passing from old electron operators ap (ag) to new 
operators Op (ap) through normal coordinates Br ■ In this way we can pass from crude 
adiabatic Hamiltonian to what is called clamped nucleus Hamiltonian and corresponding 
clamped nucleus wavefunction \1/ (r, R) . The proof that this is a canonical transformation is 
inpl. 

«P= «P + E E E Cp^Br, ...Br, ag (23) 



OO 



1 



«P = 4 + E E Zi E C'Pq-'-'^- Br, ...Br,a%, (24) 

Q k=l ri...rk 



where Br are second quantized normal coordinates. 
In short notation we can also write [|l^] 

Q k=0 

OO 

= E E C'PQ^Q 



Q fc=0 

Y^Cpqaq (25) 



a 



1 OO -| 



EE^C^T-B^a^ 



Q k=0 



{k)+ + 
PQ "Q 



Q k=0 

E C^q4 (26) 



We also perform analogous canonical transformation for phonons 

_ OO 1 

= &r + E E E <PQ '=^.1 • • • Bs.ataQ , (27) 

PQ k=0 si-Sfc 



_+ °° 1 

6. = e + EErr E CQ*5s,...Ss,4ap (28) 

PQ fc=0 si...Sfc 

The coefficients Cpg (Cpg) in eqs. (^ ^6D are determined so that ap (ap ) satisfy fermion 

anticommutation relation. The coefficients drpq ( d^pQ ) in eqs. (^ ^Sj) are determined 

- -+ _ _+ 

so that br (br ) satisfy boson commutation relation. Finally we ask fermions ap (ap ) 

_+ 

to commute with bosons br [b^ )• This means that we can write similarly as in (^ the 
total wave function \l/(r, R) as a product of fermion wave function ipk{i^, R) and boson wave 
function Xk as $(r, i?) 

^r,R)=Mr,R)Xk{R), (29) 

It is easy to show that we have two invariants of transformations, namely number operator 
of fermions 

N = N. (30) 

and normal coordinate 

B = B (31) 



The next step is that we find inverse transormations to (p5| - P8D and substitute these inverse 
transormations into eq. (p2D and then applying Wick theorem, we requantize the whole 
Hamiltonian ( [T^ ) in a new fermions and bosons [0. This leads to new V-E Hamiltonian 
(we omit - on the second quantized operators) 

H = Ha + Hp (32) 

where 
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Ha = Enn{B) - E^^Ub) - V^^\b) + Y.hRs{B)CRiC 

RSI 

+ ^ i'^RTsu ~ v%,stu)CriCsiCtjCuj 

^ RSTUIJ 

+ X! ^rsCrpCsqN 

PQRS 

+ X! i'^RTSU ~ '^^RSTu)^ RP^ SqCtiCuiN 
PQRSTUI 

+ 51 v^uvw^tpCuqCvrCwsN 
pqrstuvw 



PQr 

+ ^ (f^rA/)^ + X! ^^'^ (drPAdrQA — drPldrQl) N 

Air PQAIr 

+ Y Tl UJrdrPsdrQRN 
PQRSr 

If we introduce the following quantities 



EscF — X! ^rsCriCsi 

RSI 

+ 9 XI {^RTSU ~ ^RSTu) CriCsiCtjCuj 
^ RSTUIJ 

the new Hartree-Fock operator / with the matrix elements 



fpQ — 51 ^RsCrpCsq 

RS 

+ {^RTsu - 'Vrstu) CrpCsqCtiCui 

RTSUI 

and the new two-particle integral 



'VpQRS — Y '^tuvw^tpCuqCvrCws ■ 

TUVW 

We can rewrite our Hamiltonian Ha (^1) to the form 
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dpdQ 



ttpOjQClSOjR 



dpdQ 



a^a^asan 



Ha 



P p(2) T/(2) 



NN 

PQ 



ttpaq 



1 



'VPQRsN aj,a^asaR 



PQRS 



(38) 



Here in eq. (||) Epq IpqN 



is new Hartree-Fock operator for a new fermions (p5|), 
(^), operator | J2pqrs vpqrsN a^OqasCLR is a new fermion correlation operator and Escf 
is a new fermion Hartree-Fock energy. Our new basis set is obtained by diagonalizing the 
operator / from eq. (0)- The new Fermi vacuum is renormahzed Fermi vacuum and 

new fermions are renormahzed electrons. The diagonalization of / operator ( p6D leads to 

_+ 

coupled perturbed Hartree-Fock (CPHF) equations |]16|-[T8|]. Similarly operators hr ( ) 
corresponds to renormahzed phonons. Using the quasiparticle canonical transformations 
and the Wick theorem the V-E Hamiltonian takes the form 



H = Ha + B 



B 



(39) 



where 



and 



Ha = Emm + E 



NN "T -^SCF 

p 



CLpdp 



oo [n/2] 

n=l fc=0 

oo [n/2] 

+ E E E ^PQ~''^ N 

n=l k=0 PQ 



1 oo [n/2] 

+ ^ E E E -%^Rs'^ ■ B^"-^'^^ N Ua^qasaR 



n=0 k=0 PQRS 



+ E ^^r {drAlf + E {h^drPQ + drQpbr) N 

PQr 



Air 



apttq 



+ E {drPAdrQA — drPldrQl) N 

PQAIr 



apttQ 



+ E ^ ^rdrPsdrQRN 
PQRSr 



O'rO'qO'SO'R 



(40) 



(41) 
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As we have shown in ||T2|,|T9[ this quasiparticle transformation leads from crude adiabatic 
to adiabatic Hamiltonian. This Hamihonian (|39| ) is adiabatic Hamiltonian. Note that the 
force constant for harmonic oscilators is given as second derivative of Escf at point Ro ■ 
We shall call the corresponding phonons as adiabatic phonons. 

III. DIABATIC CANONICAL TRANSFORMATION 

In previous part we developed canonical transformation (through normal coordinates) 
by which we were able to pass from crude adiabatic to adiabatic Hamiltonian. We started 
with crude adiabatic molecular Hamiltonian on which we applied canonical transformation 
on second quantized operators 

ap=J2CpQiB)aQ (42) 
Q 

~a'p=Y.CpQ{Bya% (43) 
Q 

b=br +^drPQ{B)a~l;aQ (44) 

PQ 

K = K + Y.drQp{Bya^paQ . (45) 

PQ 

_ _+ _ _+ 

where operators ap [ap) corresponds to fermions and operators h {h^) to bosons and 
B = h + is the normal coordinate. The coefficients Cpq {drpq) can be found from the 
solution of CPHF equations. We also found that adiabatic corrections can be calculated as 
perturbation corrections, which mean that we expect that adiabatic corrections represents 
small perturbation. The situation can be more complex if we cannot treat non-adiabaticity 
as a perturbation. This is the case when non-adiabaticity can cause strong coupling between 
two or more electronic states. In order to treat such situations we can procede in a way in 
which we generalize transformations (^3|)-(^Bp. In these equations the expansion coefficients 
C and d were functions of normal coordinates B = b + b'^. 
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The generalization can be done in a way that these coefficients are some general functions 
CpQ{b, b^) and drPQ{b, b^) of b and 6+ operators. We can expect that these coefficients will 
be not only the function of normal coordinate B = b + b^ but also the function of momentum 
B = b — b^. Therefore general transformations will have the form p0| , p!9 | 



ap=Y.CpQiB,B)aQ (46) 
Q 



ap-- 



j:CpQ{B,B)^a^ (47) 
Q 



b=br + J2 drPqiB, B)a^aQ (48) 
PQ 



b, = b+ + Y,drQp{B,B)+a+paQ (49) 
PQ 

Such transformations would be rather complex, therefore we try simple approximation 

C{B,B) = C{B).C{B) (50) 

Further we can procede similarly as in the case of adiabatic approximation. We shall 
not present here the details. These are presented in 0,19]. We just mention the most 



important features of our transformation ( ^6H50| ). Firstly, when passing from crude adiabatic 
to adiabatic approximation the force constant changed from second derivative of electron- 

(2) (2) 

nuclei interaction Ug^p to second derivative of Hatree-Fock energy Eg(jp . Therefore when 
performing transformation (HUlEn) we expect change of force constant and therefore change 



of the vibrational part of Hamiltonian 

Hp = EkUB) + Ep,t{B) (51) 

The potential energy is determined by the quadratic part of the nuclear energy Ejf^{B) 
as well as by some potential energy V;^ '{B) which is a quadratic function of coordinate 
operators and has its origin in the interaction of the electrons with the vibrating nuclei. 
Therefore we have 
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EpotiB) = E^^\{B) + V^^\b) (52) 

In the case of kinetic energy term this was identical with the kinetic energy of the nuclei in the 
case of adiabatic approximation. In the case of the breakdown of adiabatic approximation 
we have to remember the finite mass of electrons and therefore to introduce more general 
kinetic energy term. Therefore, we add to the kinetic energy of the nuclei T;v(-B) some other 
yet unknown term which will be the quadratic function of momentum operator 

EuUB) = TN{B) + w'i^\B) (53) 

The total vibrational-electronic Hamiltonian 

H = Ha + Hb (54) 

will have the form 

Hs = T^^\b) + wj^\B) + eSUb) + V^^\b) (55) 

and 

Ha = Enn{B) - E^^\{B) - V^^\b) - wS\b) 

P,Q ^ PQRS 



Secondly, coefficients Cpq and Cpq are determined through equations [^| 



■"PQ + (^P ~ ^q)C'pq + X! [(^P/QA ~ ^PIAq)CaI ~ i^PAQI ^ ^PAIq)CiA 

AI 



huJrCpQ = EpSpQ (57) 

(^P ~ ^q)^PQ + X! [(^P/QA ~ ^PIAq)^AI ~ i^PAQI ~ ^PAIq)^IA 

flUJrCpQ = EpdpQ (58) 



AI 
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where hur is the new non-adiabatic phonon given by 

1. 



(J 



The expressions for extra terms vI^\b) and W^'{B) in (|55D are given as 



(2), 



yrs 



lA 



and 



AI 



This means that the resulting vibrational frequency u depends explicitly on coefficients C 
and CpQ. 

Finally fermion part of Hamiltonian will be given as 



Hp = Hp + i^p' + Hrp" + Hr 



For the ground state energy we get 



Hp = + EgcF + E ^^tUCaI^ - \C\i 



r \2\ 



Air 



One-fermion part will be 



Hp' = ^6pN[aj,ap] + E 

P PQr 



. A 



KypA^QA ~ ^pa^qa)~ 



PQr 



+ E [(hiP) -p{PWp + ih{Q) ~ p{Q))e^Q*] C},QN[ataQ] 



PQr 



- E [i'^PiQA - Vpiaq)Ci*a + (vpAQi - Vpaiq)Ca*i] Nla^aq] 

PQAIr 



Two-fermion part will be 
1 



Hp" = 7^ E VpQRsN[(^p(^Q(^saR] 



PQRS 
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+ Y: ^^r{C^PRCsQ-C}>RC7Q)N[ata^asan] 

PQRSr 

-2 £pCs*QN[a^a-^asaR] 

PQRSr 



+ 2 E E[' 

PQRSTr I / 



"^PQTS^TI "^PQTI^TS + i'^TQSI '^TQISJ'^PT 



^RI 



'^TIRS^QT '^TQRS^IT + ("^/QTS '^^QSTJ^TR 



^RA 



(65) 



- E \^tarsCqt - vtqrsC'j^t + i^AQTS " i'aqst)Ctr C'^/p \ N[apa'^asaR\ 
A ) 

Three-fermion part will be (as a result of transformation (^) the three fermion term appears) 
Hf'" = -2 Y (^pqvtC'rs - VvQSTCpv)CuRN[aj>a:l^a^auaTas] (66) 

PQRSTUVr 

The bosonic part of Hamiltonian Hp is not given in a diagonal form. To bring it to diagonal 



form as in eq. (|59D we can proceede as follows. 



Hb — Tn + e'^nn 



/ . EscF ~ E ^11 +9 + ^11 

r,s / L ^ \ 



+ E (//?/ Sri + h LJr Crj] C ri 

RI ^ 

/A'' ^ A'' ~s \ 

+ [fm Sri +hu;s Crj) Crj BrBs 



AI 



(67) 



Our aim is to bring this Hamiltonian into diagonal form. We can extract adiabatic part 
{huj'^) and we get 



where 



^;. = 2E E 

r,s I RJ I 



UrI " )i.CRj — C^/) + ^ ^rCjijCRj 



BrBg 



(68) 



(69) 
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(70) 



AI 



where Cj^j is identical with Cri coefficients from adiabatic transformation eqs. ( p3| , p^). If 
we substitute in eq. ( |69D and (|70]) for Br = + br and Br = K — b^ we get for ( p8| ) the 
expression 

1 



where 



Arsb^bs + ^BrMbj + brbs 



(71) 



Brs = '^{F^s + F, 



(72) 



pi _|_ pi _ P2 _ p2 

rs ' sr rs sr 



(73) 



Diagonahzing the above Hamiltonian we obtain diabatic frequencies 



1, 



Hb = Y1 h^tiPtbr + -) 



(74) 



Hamiltonian ( [7T| ) has a form of quadratic Hamiltonian [^,0 and can be diagonalized by 
Bogoljubov transformation, which leads to the condition 



det 



A-u^l B 
-B -A - u'^i 







(75) 



Secular equation ( [75[ ) gives us diabatic phonons hUr- 
If we look at eq. ( pj) we see that we have corrections due to non-adiabaticity to one-particle 
part as well as to two-particle part. We see the hierarchical structure of our Hamiltonian. If 
the non-adiabatic coupling is small i.e. C goes to zero and we have adiabatic Hamiltonian. 
If this coupling is strong we cannot use adiabatic approximation but we have to work with 
full V-E Hamiltonian (M). 
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IV. CALCULATIONS 



In order to compare our approach with other approaches deahng with adiabatic correc- 
tions we perform simple model calculations for adiabatic corrections to ground state energy. 
We start with adiabatic Hamiltonian (^). We now perform the following approximation. 
We limit ourselves to finite orders of Taylor expansion of the operators H'^ and H'^ . We shall 
use similar approximation as in |^ . The diagrammatic representation of our approximate 
Hamiltonian will be 



H = Ho + H = E%j^ + E%f;p epN 



dpCip 



+ J2 hcOr[btK + ^) (76) 



Here should be included image. gif- 



The adiabatic corrections to the ground state of H2, HD, and D2 we shall calculate using 
second-order Rayleigh-Schrodinger many-body perturbation theory (RS-MBPT) and our 
Hamiltonian (|76D . If we assume that we know the solution of the unperturbed Schrodinger 
equation 
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-f^o I V^o > = Eo\ipo> , (77) 
where Hq is the unperturbed Hamiltonian + where 

HI = E"^j^ + E°scF +Y.^P^[^>p] (78) 

p 

and 

H'^B = [KK + ^) (79) 

The perturbed (exact) Schrodinger equation will read 

H \ > = £\-^ > , (80) 



where H will be our Hamiltonian (ffq). The perturbed energy E will be given through the 
RS-MBPT expansion as 

£ = Eo+ < (po\ H' \ ifo > + < (po\ H'QoH' \(po > + ■■■ , (81) 
where H' is the perturbation and Qo is the resolvent 

= E^^^ (82) 

Since our sets of boson creation and annihilation operators and fermion creation and an- 
nihilation operators commute we can write our unperturbed wavefuntion \ > as the 
product of the fermion state vector | ipo > and the boson state vector \ Xo > , i-e. 

I ifo >=\ iIjo> \Xo> (83) 

Further we want to study the nonadiabatic corrections to the ground state. Therefore 
I tpo > will be the unperturbed ground state wave function (we shall use Hartree-Fock 
ground state Slater determinant -Fermi vacuum) and | Xo > will be boson ground state 
-boson vacuum | >. 

I Xo >= I > . (84) 
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The exact ground state energy will be given by perturbation expansion (up to the second 
order) 

£ = < i>o \ < \ Ho\0 > \ ^/'o > 
+ < ^0 I < I H' \ >| V^o > 

+ < 7/^0 I < I H'QoH' I >| ^0 > + • • • (85) 
Substituting for Ho = + H% from (^, |7^) into the first term in (^5]) we get 

< Vo |< I //o I >| ^0 > = ^^7V + E°scF + E \^^r (86) 
We shall not present all terms for perturbation corrections from the right hand side of 



eq. dSSf ). There are corrections which corresponds to electron correlation, anharmonicity 



corrections and adiabatic corrections [^. We shall pay attention only to adiabatic correc- 



tions given through second-order term in eq. (|85|) . Analyzing diagrammatic contributions 
through the Hamiltonian ( [76D we find that the adiabatic corrections are given through the 
second up to fourth term in eq. (|4l|). From these terms we calculate only contributions from 
the first and the second term which are given through the first order of Taylor expansion and 
these terms are used in second order RS-MBPT. We obtain the following simple expressions 

< I < I ij' I > I -i/^o > ~ 

Y ^^r [CaI + < ^(0) \I> 
Air ^ 

+ E ^'^r {Cai + < A(0) |f >)' . {si -SA- TiUrY^ 



Air ^ 

f l\r ST \ 

+ E ^^^r \CaI^< ^(0) 1/ > ) . [ei-EA^ hUr)'^ 
Air ^ 

^Y.'^^r [CaI + < ^(0) 1/ > 
Air ^ ' 

^2Yn^ujl{cAi^<m\l>] .{ei-EAY^ (87) 



Air 
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For the notation see [^,^. We believe these three terms on right hand side of (|87|) rep- 
resents the dominant contributions to adiabatic corrections. The last formula is valid due 
to the eq. (|^). This formula was used to calculate the adiabatic corrections to the ground 
state energy of the if2, -D2, and HD, molecules. Results in cm~^ are presented in Table I. 
We see that the structure of this formula is similar to eq. (|^ except that we use RS-MBPT. 
We also used the same approach to calculate the adiabatic corrections to the energies of 
the first vibrational transitions for the same molecules Analyzing eq. (^) we can see 
that the first term on the right hand side of eq.(^) is always positive and correspond to 
largest contribution, while the second term on right hand side of eq. (^) is always negative 
and represents smaller contribution than the previous term. Therefore we can expect that 



the eq. (|87D should converge to the true value of adiabatic correction from above. This 
also explain the larger values for H2-, D2, and HD adiabatic correction obtained through 



eq. ([87| ) than true value obtained by Wolniewicz p^. Another source of difference can be 
basis set used and also the contribution from other terms in (Pl|). In our calculation we 
have obtained for H2 value of 136.89 cm~^ using Roos Augmented Triple Zeta ANO 
basis set. Handy using basis set of similar quality as used in our calculations obtained 
value 101 cm~^. Similar value of adiabatic correction was obtained also in an older study by 
Wolfsberg Kutzelnigg in his paper using wave function expansion consisting 1200 
functions obtained Wolniewicz p4| value 114.591 cm~^. It is apparent that calculations of 
adiabatic corrections are strongly basis set dependent. 



V. CONCLUSIONS 

In this article we performed simple model calculations of adiabatic corrections for ground 
state energy of the if 2, HD, and D2 molecules. The corrections were derived through canoni- 
cal transformation applied to crude adiabatic molecular Hamiltonian. These transformations 
mix together electrons and phonons (normal coordinate canonical transformation) leading 
to adiabatic molecular Hamiltonian. Using second quantization formalism and many-body 
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diagramatic perturbation theory and splitting the adiabatic Hamiltonian into unperturbed 
part and perturbation we derived the formulae for adiabatic corrections. The results were 
compared with the obtained by different approaches by Wolniewicz [E4| and recently by 



Kutzelnigg ^7j. The quasiparticle canonical transformations were then generalized in a 
way that electrons and phonons are mixed not only through the normal coordinate but 
also through the momenta. This canonical transformation leads to non-adiabatic molecular 
Hamiltonian (motion of electrons does not follow the motion of nuclei, the electrons are 
phase shifted with respect to nuclei). One can clearly see that the electronic and vibrational 
motion cannot be separated. The mixed system behaves as one whole quasiparticle (mixed 
electrons and phonons through the last (momentum) transformation behaves as a renormal- 
ized fermions and a different mixture of electrons and phonons leads to renormalized bosons). 
This in some extent analogous to the introduction of quasiparticles in the solid state theory, 
where the "bare" electron interacting with quantized lattice vibrations is renormalized to 
"absorb" some part of this interaction, and this quasiparticle is known as a polaron. We 
were able to derive equations for non-adiabatic C coefficients, which permits us to calculate 
the so called mass polarization terms and thus non-adiabatic phonons. It is interesting that 
the quasiparticles preserves some interesting features known from pure electronic molecular 
Hamiltonian calculations e.g. we can speak about orbital energies of a new quasipartcles, 
correlation energies of of a new quasipartcles, Hartree-Fock energy (holes and particles), 
etc. Further very important property which follow from the last canonical (momentum) 
transformation is that we clearly see that in the case of electronic quasidegeneracy when C 
coefficients are non-negligible (non-adiabaticity is not a small correction) we should work 
with a full non-adiabatic Hamiltonian. In the case when non-adiabaticity is a small cor- 
rection and C coefficients are negligible, we can work with adiabatic Hamiltonian (we have 
only C coefficients through CPHF equations) and only if system is perfectly separable that 
even C coefficients are negligible we can work with purely electronic Hamiltonian. 
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TABLES 



TABLE L Adiabatic corrections (in cm ^) for H2, D2 and HD molecules. 





Wolfsberg |26| 


Wolniewicz ||24| 


Kutzelnigg ||27| 


Handy llO| 


Present method, eq.(|87|) 


H2 


101.3 


114.591 


114.591 


101 


136.89 




50.7 


57.296* 


57.296* 


50.5* 


68.62 


HD 


76.0 


85.943* 


85.943* 


75.8 * 


102.67* 



*These adiabatic corrections were obtained by reduced mass rescaling from H2 corrections 
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